import numpy as np
import math
import matplotlib.pyplot as plt

#omiga = 9
#file_name1 = "C/data/Excitation_C_2s22p23P-2s22p3p3D_ref.dat"
#file_name2 = "C/data/Excitation_C_2s22p23P-2s22p3p3D_lxcat.dat"

#file_name1 = "C/data/Excitation_C_2s22p23P-2s22p21D_ref.dat"
#file_name2 = "C/data/Excitation_C_2s22p23P-2s22p21D_lxcat.dat"

#file_name1 = "C/data/Excitation_C_2s22p23P-2s22p3s3P_ref.dat"
#file_name2 = "C/data/Excitation_C_2s22p23P-2s22p3s3P_lxcat.dat"

#file_name1 = "C/data/Ionization_C_to_C+1.dat"
#file_name2 = "C/data/Ionization_C_to_C+1_lxcat.dat"

#file_name1 = "C/data/Excitation_C_2s22p23P-2s2p35S_ref.dat"
#file_name2 = "C/data/Excitation_C_2s22p23P-2s2p35S_lxcat.dat"

#file_name1 = "C/data/Excitation_C_2s22p23P-2s22p21S_ref.dat"
#file_name2 = "C/data/Excitation_C_2s22p23P-2s22p21S_lxcat.dat"

#file_name1 = "C/data/Excitation_C_total.dat"
#file_name2 = "C/data/Excitation_C_total.dat"

file_name1 = "C/data/Ionization_C_to_C+1_ref.dat"
file_name2 = "C/data/Ionization_C_to_C+1_lxcat.dat"


cs1 = np.loadtxt(file_name1)
cs2 = np.loadtxt(file_name2)

print("cs1 shape: ", cs1.shape)
print("cs2 shape: ", cs2.shape)

plt.plot(cs1[:,0], cs1[:,1], label = "1")
plt.plot(cs2[:,0], cs2[:,1], label = "lxcat")

plt.legend()

plt.savefig("compare_cross_section.svg")
plt.show()
